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ABSTRACT 

Aims. We investigate the dissipation of the bulk kinetic energy of a relativistic jet at different distances from the central power-house and 
analyse in detail how the dissipated energy is radiated away. 

Methods. We assume that the location of the dissipation region is a function of the bulk Lorentz factor T of the jet, being closer to the center 
for smaller T. This assumption is naturally fulfilled in the internal shock scenario. The dissipated energy is partially used to accelerate electrons 
and to amplify the magnetic field. This process creates a source inside the jet (blob). Such blobs may efficiently produce synchrotron and 
inverse Compton emission. 

Results. We find that even if the blobs or shells responsible for the blazar activity carry the same energy (in bulk kinetic form), the fact that 
they move at different T can produce dramatic variations in different bands, even if the bolometric luminosity is instead very similar. This is due 
to the relative importance of the synchrotron, self-Compton and external Compton radiation processes, which greatly changes by changing T 
and the compactness of the source, even if the total radiated energy is constant. We then find that the jet can produce most of its radiative output 
at small distances from the putative black-hole and its accretion disk, if this implies a low level of emitted MeV-GeV flux. Our findings, which 
we apply for illustrative purposes to the blazar 3C 454.3, will be easily testable by the coming y-ray satellite, such as AGILE and GLAST. 
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1. Introduction 

Our knowledge of the spectral energy distribution (SED) of 
blazars, at MeV-GeV energies, has been built upon the results 
of the EGRET instrument onboard the Compton Gamma Ray 
Observatory, which shows that blazars as a class are powerful 
y-ray emitters, and that the entire SED is characterized by two 
broad peaks, with the high energy one becoming more domi- 
nant for more powerful blazars (Fossati et al. 1998 1. This al- 
lows to define a blazar sequence, whose main parameter is the 
bolometric observed luminosity, dictating the overall spectral 
appearance of the SED (Ghisellini et al.[TP98 1: in more power- 
ful objects, in which the radiative cooling is stronger, the emit- 
ting electrons have relatively small typical energies, explaining 
why in these sources the peak frequencies of the two peaks are 
smaller than in low power blazars. Since in these sources the 
stronger cooling is mainly due to the presence of external radi- 
ation (e.g. the broad emission lines), this scheme also explains 
why in powerful sources the high energy peak dominates. This 
blazar sequence has been obtained by averaging the data of a 
relatively large number of blazars (more than one hundred), and 
not all of them have been detected by EGRET. It then accounts 
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for the fact that a sizeable fraction of blazars went undetected 
by EGRET. On the other hand, this average does not entirely 
account for the fact that one specific source can spend most 
of the time at MeV-GeV fluxes which are below the detection 
threshold. Indeed, only a fraction (roughly 1/4) of bright radio 
blazars were detected by EGRET. This can bias the representa- 
tion of the average SED in the sense that the high energy flux 
might better represents the high state of the source, rather than 
its time-average SED. 

A related issue concerns the possibility of independent vari- 
ability of the two peaks. Judging from the (few) examples we 
have, we know that variability at high energies (fluxes belong- 
ing to the high energy peak) is accompanied with flux variabil- 
ity at smaller (i.e. UV, optical, IR, but not radio) frequencies 
(with the rare but interesting exceptions of the so-called "or- 
phan flares", i.e. TeV flux variations not accompanied by si- 
multaneous X-ray variations; Krawczynski et al. 20(3). 

We then wonder if it is possible that a blazar flares in its 
synchrotron part of the spectrum, but not in the MeV-GeV 
band. As an illustrative example, consider the blazar 3C 454.3, 
detected by EGRET at the beginning of nineties (Hartman 
et al. I1999> . when the source was rather faint in the opti- 
cal, and by BeppoSAX (in 2000, Tavecchio et al. l^UOll . The 
BeppoSAX spectrum, up to ~ 100 keV, was relatively faint but 
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very hard. Its extrapolation to the EGRET band was roughly 
consistent with the flux detected by EGRET in a different pe- 
riod of time. Recently, 3C 454.3 showed a huge flare in the 
optical (Fuhrmann et al. 2006 Pian et al. 2006 1 and in the 
X-ray bands, up to 100 keV, as detected by INTEGRAL and 
SWIFT (Pian et al. l2006l Giommi et al. l2U06l . However, we 
have no information about the emission level in the MeV-GeV 
range during above mentioned activity. The SED of 3C 454.3 
was interpreted by Ghisellini et al. (Q]598 1 and Tavecchio et al. 
(2002 1 as due to a synchrotron and Inverse Compton model, 
produced by one-zone region of the relativistic jet, in analogy 
with other powerful blazars. The high energy peak was the re- 
sult of the inverse Compton scattering of relativistic electrons 
off seed photons coming from the broad line region (BLR), 
with the synchrotron self-Compton process playing a negligi- 
ble role. If we apply the very same model to the optical/X-ray 
flaring state observed recently, then we are led to conclude that 
also the MeV-GeV emission should show a flare as well, of 
similar amplitude to what seen in the optical . We are then 
led to conclude that the jet has indeed changed its total power, 
by a large factor. Unfortunately, this scenario cannot be tested 
immediately due to lack of the MeV-GeV observations. There 
are, however, more "economic" solutions, in which the total jet 
power changes much less, or it is even constant. One solution 
was proposed in Pian et al. (2006 1: in this paper it was pro- 
posed that the emission site, responsible for the flare observed 
in 2005, was outside the BLR, where the number of external 
photons was negligible. In this case the external Compton com- 
ponent becomes negligible, and the power is mainly emitted by 
the synchrotron and self-Compton process. With roughly the 
same jet power, we can account both for the "high EGRET" 
state (with a large dominance of the high energy peak over the 
synchrotron one) and the "high synchrotron" state, where the 
power in the synchrotron and synchrotron-self Compton com- 
ponents are roughly equal, and equal to the previously observed 
"high EGRET" state. 

In this paper we explore yet another possibility, in which 
the emission site is within the BLR, but the importance of the 
external Compton component is much reduced because of a 
reduced bulk Lorentz factor F. In fact, the importance of this 
component depends on the radiation energy density U' ext of 
the seed external photons as seen in the comoving frame, and 
U' ext oc r 2 . It is then possible, in principle, to reduce the external 
Compton component and at the same time to enhance both the 
synchrotron and the self-Compton ones. To be more quantita- 
tive, we have also to specify not only the F-factor, but also the 
size of the region, the magnetic field value, the density of the 
particles and their energy spectrum. There is a specific scenario 
where all these quantities are all dependent in a specific way to 
the value of F, and this is the internal shock model. We will dis- 
cuss it in detail in the next section, but the basic idea is that in 
this scenario the central engine is working intermittently, pro- 
ducing shells of matter moving at different speeds (and hence 
different Lorentz factors). At the start they are separated by 
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same distance AR, and there is always some probability that 
a later shell is faster then the previous one. The fastest shell 
will then catch up the slower one, at a distance of the order of 
^coii ~ r 2 AR. Therefore, if, on average, the two shells are pro- 
duced with T factors that are different, but relatively small, R co \\ 
will be correspondingly small and the emission region will be 
more compact, and the particle density and the magnetic field 
will be greater. If, in addition, we also assume that the kinetic 
and magnetic energy carried by the shells is always the same 
and that the dissipated energy (in the collision) is also the same, 
we have that all quantities depends on F only. 

The internal shock scenario then offers a nice possibility to 
explore and quantify the consequences of the idea that a typi- 
cal relativistic jet, even if on average produces blobs of equal 
energies and therefore works, on average, at a constant power, 
yet it can produce spectacular flares in different energy bands. 
However, we would like to stress that the main idea explored 
in this paper is not to test the internal shock scenario, but to in- 
vestigate what happens if the jet, carrying on average the same 
amount of energy, dissipates it at different distances from the 
black hole, characterized by different T-factors (smaller T at 
shorter distances). Our study also allows us to answer a more 
general question: Is it mandatory that the jet always produces 
most of the radiation we see at relatively large distances from 
the black hole? If the jet dissipation always results in the pro- 
duction of a dominant MeV-GeV component the answer is yes. 
As discussed in Ghisellini & Madau ( 1 996 1, if high energy pho- 
tons are produced too close to the accretion disk and its X-ray 
corona, they would be absorbed by photon-photon collisions, 
producing electron-positron pairs. These pairs would Compton 
scatter the UV-optical radiation from the disk, reprocessing 
the energy originally in the MeV-GeV band into the X-ray 
band. The observed "valley" between the two broad peaks of 
the blazar SED would then be filled by this reprocessed radi- 
ation. But what happens if the jet, close to the disk does not 
produce a too prominent MeV-GeV component? 



2. The model 

In our model we assume a relativistic jet that contains electrons, 
protons and tangled magnetic field. Different sub-structures 
of this jet (hereafter shells) may travel with different veloc- 
ities. Therefore, some of the shells may collide and gener- 
ate relativistic shock waves. The shock, through the first or- 
der Fermi acceleration, may increase significantly the energy 
of the electrons. In other words, a fraction of the bulk kinetic 
energy of the colliding shells is used to accelerate electrons 
up to highly relativistic energies. The relativistic electrons es- 
cape from the shock front into the downstream region of the 
shock, where they lose their energy through the synchrotron 
and the inverse-Compton processes. These are the main as- 
sumptions of the so-called internal shock scenario, proposed 
for the first time by Rees (1978 1 and applied successfully for 
gamma-ray bursts (see e.g. Meszaros 2006 for a recent review) 
and blazars (Sikora et al. l 19941 Begelman et al. 1994 Ghisellini 
IT9991 Spada et al. l200T1 Guetta et al. 2004). 
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2.1. Set up of the model 

Constructing our model we make two important assumptions. 
The first concerns the efficiency of the central engine that forms 
and accelerates the shells. We simply assume that the efficiency 
of the central engine is constant. This assumption is not a key 
ingredient of our model, but it allows to see if even if the jet 
is working with the same efficiency, we can account for large 
variation of the flux in different bands, as observed. According 
to this first assumption, the bulk Lorentz factor T — 1 / yjl - 1 
of the shell is inversely proportional to the total mass of the 
shell, so that the total energy given to the shell by the central 
engine is constant: 



MF = const. 



(1) 



where M N p m p (N p is total number of protons of mass m p ). 

Consider now a collision of two shells (a and b) of bulk 
kinetic energy T a M a and V b M b respectively. The collision must 
satisfy energy 



Y a M a + TbMb = r s (M a + M b + E'/c 2 ), 



and momentum 



T a [3 a M a + T b p b M b = T s p s {M a + M b + E'/c 2 ), 



(2) 



(3) 



conservation, where T s , fi s are describing the "merged" shell 
(s) created from the colliding shells and E' is the energy dis- 
sipated during the collision, measured in the comoving frame 
of the new shell. As mentioned above, this energy can be used 
to accelerate the electrons, and it can also be used for heating 
the protons or to amplify the magnetic field intensity. Since it is 
the new shell which contains the relativistic emitting particles, 
hereafter we will call it simply "source" or "blob". The effi- 
ciency describing how much of the kinetic energy of the two 
colliding shells is dissipated can be defined as 



(4) 



T a M a + T h M b 



where E is the dissipated energy in the observer's frame. Using 
the conservation laws (Eq.[2]and Eq.[3} we can derive 



n = i- 



r s (i + a M ) 

r a (l + a T a M y 



(5) 



where % = M b /M a and ffp = T b /F a . Since our first assump- 
tion (shells of constant bulk kinetic energy) implies ar = 1/<*m> 
the above formula simplifies into 



>1 



1 



1 + % 
V2 + 2a, 2 „ 



(see e.g. Lazzati et al. 1999 1. This simple relation shows that r\ 
varies between and 22% if the contrast between the Lorentz 
factors (ap) changes from 1 to 10. 

Our second assumption concerns the amount of the energy 
transferred to the electrons during the collision. We simply as- 
sume that the electrons receive always the same amount of en- 
ergy (E e - const.). Since our first assumption implies constant 
energy of the colliding shells, the present assumption requires 



only that the contrast between the Lorentz factors of the collid- 
ing shells is always the same (ffp = const.). We are aware that 
this may be a great simplification of what we can expect in real- 
ity. However, we will see that this simple assumption is enough 
to provide a reasonable explanation for the some properties of 
blazar variability. 

The total energy transferred to the electrons is given by 



E e - YE' e — TN e < y > m e c - const., 
where N e is the total number of electrons and 



< y >= 



£7 yQ(y)dy 

r- Q(y)dy 



(7) 



(8) 



is the average electron energy, where 7=1/ -yl - ft 2 , is the 
electron random Lorentz factor and Q(y) describes the injec- 
tion of the particles in the energy range constrained by y m i n 
and y max . In the model we describe the particle acceleration as 
a continuous injection of relativistic electrons. This injection 
is simulating the escape of the particles from the shock region 
into the downstream region of the shock that forms the source. 
First, we will discuss the very simple case of a monoenergetic 
injection, later we will discuss more complex, but more realis- 
tic, injections of power law energy distributions. The duration 
of the injection is assumed to be equal to the source light cross- 
ing time (R/c). During the process the electrons that are car- 
rying the total energy described by Eq. are injected into the 
source. The evolution of the particle energy distribution N(y, t) 
inside the source is calculated from the kinetic equation: 



dN(y, t) 



dt 



C 



dy 

4 cr T c 
3 m e c 



— {Cy 2 N{y,t)) + Q{y), 
2 (U' B + U' syu + U' ext ), 



(9) 
(10) 



where U' B , U' sya , U' ext are the magnetic, synchrotron and exter- 
nal radiation field energy densities, respectively. The solutions 
of this equation for the different types of injections, that we use 
in our modelling, were first derived by Kardashev dl962> . 

For the sake of simplicity, the source is approximated as a 
spherical homogeneous blob, of radius R assumed to be pro- 
portional to the distance D to the center of the active nuclei 
(R = ifrD, where \]/ is the jet aperture angle). We also assume 
that the D depends on the bulk Lorentz factor T as D oc F 2 . This 
assumption is supported by the internal shock model where two 
shells (Ft, > F fl ) will collide when the second source (r^) will 
travel the distance 



(6) £> = 2- 



D r 2 „ 



(ii) 



where Do is the initial separation of the sources (see e.g. 
Lazzati et al. 119991 1. In other words, relatively slow shells 
(r «: 10) may collide closer to the center than faster shells 
(r > 10) that can collide only at larger distances. Note that this 
formula is valid only for constantly moving shells (no acceler- 
ation or deceleration). Moreover, this formula shows that our 
assumption (D oc F 2 ) becomes inaccurate for small F-contrasts 
(ar < 2). Finally, our assumption is supported by the fact that 
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the central engine is probably not able to accelerate almost in- 
stantaneously the jet components up to very high values of T. 
If the acceleration process requires a time comparable to the 
travel time of the central engine (e.g. 0.1 pc) then we should 
expect fast shells farther from the center than slow shells. 

To complete the general description of our model we have 
to define the evolution of the magnetic field. By analogy to the 
definition of the energy transferred to the electrons (Eq.[7), we 
assume that the total energy accumulated after the collision in 
the magnetic field, measured in the observer's frame, is con- 
stant 



E B = YE' = TU'V = rB' 2 V'/(8n) = const., 



(12) 



where B' is the magnetic field intensity and V is the source 
volume as measured in the comoving frame. In other words, the 
fraction of the total dissipated energy used for the amplification 
of the magnetic field is constant. 

Finally, some part of the dissipated energy can be trans- 
ferred to the protons. However, we neglect this process in our 
computations since protons are not assumed to radiate and 
since their energy changes slightly for different Lorentz fac- 
tors. 

2.2. Comparison of different sources 

Our main aim is to investigate how the energy transferred 
to the electrons (that is constrained to be exactly the same 
for every source) is radiated by different sources. By "differ- 
ent sources" we here mean different states of the same jet. 
Therefore, we simply compare the emission of two sources 
with different Lorentz factors (r 2 < H). The difference in the 
Lorentz factors means that the sources were created by colli- 
sions of a pair of shells with different velocities (e.g. F/, = 20 
andr fl = 10 -» Ti 13.3 andT c = 8 and T d = 4 -> Y 2 =* 5.3). 
Note that according to our assumption the contrast between the 
Lorentz factors of the colliding shells (o^ = Fft/F fl = r c /r d ) is 
constant. Moreover, we assume that the shells before the colli- 
sion contain only "cold" particles that cannot emit. Therefore, 
it is not mandatory to describe the physical parameters of both 
sources: it is sufficient to constrain the physical parameters of 
the first source using observations of the first state and derive 
the parameters of the second source (second state) choosing T 2 
and using our assumptions 



,B n = B, 



•3.5 



, < 7 >2 



N e ,l < J >1 



N, 



■(13) 



e.2 



Finally, to complete the description of the sources it is neces- 
sary to describe the injection function. 

2.3. Monoenergetic injection 

In the first approach we assume continuous monoenergetic in- 
jection of relativistic electrons 



G(r) = Qi Dirac(y - y t ), 



(14) 



where Qi is the number of injected particles per time unit and 
ji describes the energy of the injected electrons. This approach 



helps to simplify the description of the model (e.g. < y >= y,) 
and allows to derive simple formulas that describe the evolu- 
tion of the observed emission. According to the solution of the 
kinetic equation, the N(y) distribution is a power law 



N(y) 



C 



y " for y min < y < y t ; y 

IT 



Ji 

\ + c 7i r 



(15) 



where n — 2 in this particular case, T — R/c is the duration 
time of the injection assumed to be equal to the source light 
crossing time. The total number of injected electrons is 



N. 



-vrj 



Q(j)dj, 



(16) 



that gives N e = VT Qi in the particular case of a monoen- 
ergetic injection. We assume similar number of electrons and 
protons inside the shells before the collision (e.g. N e « - N PA + 
N e ,b — Npj,). However, only some fraction of the electrons 
may be accelerated. Moreover, we assume that the total num- 
ber of accelerated electrons (total number of the injected elec- 
trons) is proportional to the total number of protons inside the 
source. According to our main assumptions (T a M a = TbMb = 
r c M c = TjMj — const, and T\E'\ = T 2 E' 2 = const.) we obtain 
Y\{M a + M b ) - T 2 (M C + M d ). If we neglect the contribution of 
the electrons to the total mass of the shells then the total num- 
ber of protons and therefore also the total number of electrons 
in different sources scales as 

N 2 =N l -^-. 

1 2 

This gives 



(17) 



Qui = Qi, 



V' 2 T' 2 T 2 UuX \Y x 



(18) 



where the injection time (T') and the source volume (V) scale 
as 



(19) 



Moreover, since y t =<y> the transformation of the total num- 
ber of the particles gives y- u2 = yi,{YilT\. The electron density 
inside the source after the end of the monoenergetic injection 
is defined by 



K 



Qi 3m e c 2 



Q, 



C 4a T c U'+U'+U' 



(20) 



where the magnetic field energy density (U' B ) determines the 
synchrotron cooling during the injection, the synchrotron ra- 
diation field energy density (U' syn ) determines the cooling due 
to synchrotron self-Compton emission (this is correct only in 
the Thomson limit) and the radiation field energy density of the 
external medium, measured in the comoving frame of a source 
(U' en ) determines the cooling due to inverse-Compton scatter- 
ing of external photons. We can relatively easily describe the 
change of U' B and U' ext as a function of F. However, the pre- 
cise description U' is quite complex and requires a numerical 
approach. Therefore, for deriving some analytical formulas de- 
scribing the emission processes, we will focus on two simple 
and opposite cases. 
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First, we assume the dominance of the EIC scattering 
(U' ext » U' B and U' ext » U'^). Moreover, U' ext is assumed to 
be constant in the observer's frame with a good approximation 
if this external radiation is produced by the Broad Line Region 
(BLR) (see e.g. Ghisellini & Madau 1996 1. Therefore, the ex- 
ternal radiation field density in the comoving frame for sources 
with different Lorentz factors scales as 



^ext,2 



If 

u ext,l 



r 2 



(21) 



This gives the following transformation of the particle density 
K 2 = K, 1^1 , (22) 

Knowing how all physical parameters of the source (/?', B', 
K, T, U' ext ) scale with F, we can derive the approximated for- 
mulae describing the different types of emission, and how these 
depend on the bulk Lorentz factor. 

First, consider the synchrotron emission where the emissiv- 
ity at given frequency depends on the particle density and the 
magnetic field intensity as 



(23) 



where a = (n - l)/2 is the energy spectral index. In the case 
of mono-energetic injection, a = 0.5. The observed thin syn- 
chrotron flux is given by 



F syn (v) oc 6 



3+ff D 2 



(24) 



where 6 = [T(l —BcosO)] 1 is the source Doppler factor, 9 is 
the viewing angle and I^ yn (v') K R 'fsyJ v ') K R'KB n+a v'- a is 
the intensity of the thin synchrotron emission as measured in 
the comoving frame. For simplicity we here neglect the self 
absorption regime (we will mainly investigate the synchrotron 
emission in the optical and X-ray energy bands where self- 
absorption is negligible). Moreover, in order to simplify the 
description even more we assume F = 6. This is strictly true 
only if the viewing angle 9 = 1/F [rad]. This assumption en- 
ables possibility to describe the flux as a function of the Lorentz 
factor only. Strictly speaking this implies that when we com- 
pare the fluxes produced by sources with different Lorentz fac- 
tors we compare these sources in different astrophysical ob- 
jects with different viewing angles. However, this approxima- 
tion does not change qualitatively our results in comparison 
to the case where we fix the viewing angle (this is mandatory 
when we study the same astrophysical object) and modify the 
value of the Doppler factor. Note that for 9 = the Doppler 
factor 6 ~ 2F. We use this approach only for the tests presented 
in this subsection for illustrating the relations between differ- 
ent types of the source emission in a simpler way. Using the 
above assumptions we can write F syn (v) oc r 35 R 3 KB' l5 v~ 05 
that gives the following scaling 



^.syn,2(v) = /\syn,l(v)| — 



-6.75 



(25) 



This shows that the synchrotron specific flux increases when 
T decreases. This is because the increase of the magnetic field 



intensity and the particle density can easily compensate the de- 
crease of source volume and the decreased Doppler boosting. 

In a similar way we can describe the evolution of the self- 
Compton emission where the emissivity depends on the par- 
ticle density and the intensity of the radiation field j' ssc (v) oc 



Ki: yn (v , y. 



issc(v') <* KI'W) <* R'K l B 



,/p-2 D /1.5 /-0.5 



2d/1.5„-Q.5 



The observed flux is F ssc (v) oc r"fl 4 # z fi 
cording to our previous assumptions scales as 



(26) 



which ac- 



(27) 



Note the stronger dependence on T with respect to the syn- 
chrotron flux. 

Finally, the external inverse-Compton (EIC) flux is given 
by Feic(v) = f E ic,ov~ - 5 oc r 3 - 5 R 3 KF ext V-°- 5 . The (frequency 
integrated) intensity of the external radiation field in the co- 
moving frame (7 ext ) is proportional to F 1 5 . To understand this, 
assume a source with a constant size, and particle density, sur- 
rounded by a constant external radiation field. Change only F. 
The total flux (F{™r.) produced by such a source is proportional 
to F 2 <!> 4 , and, according to our assumption (r = 6), F^L oc T 6 . 
Assuming a power law spectrum of the emission we can write 



^EIC 



Since 



I 



r -0.5 j 

Feic,ov dv 



^EIC,0V max . 



(28) 



oc SF ~ F 2 , this relation gives Feico k 
Moreover, for the source where R and K do not depend on V 



we have Feic(^) k T IL 



oc r that gives I' ext (v') oc r 



Using this relation and the previous definitions we obtain 

Feic,2(v) = F mc ,i(v), (29) 

as long as we compare the fluxes at a given frequency in 
the power law part of the spectrum. This means that our two 
sources with different Lorentz factors have the same specific 
flux if the total energy transferred to the particles is the same 
for each source and the physical parameters of the sources are 
scaling according to our prescriptions. This result is valid only 
for the sources where the emission is dominated by the external 
inverse-Compton scattering. 

Consider now the opposite case where the cooling is dom- 
inated by the synchrotron emission (U' B » U' ext and U' B » 
U' syn ). The magnetic field energy density transforms as 



that according to Eg. 1201 gives 



(30) 



(31) 



This scaling differs significantly from the result obtained in the 
previous case (Eg. I22>. where U' ext was dominant. Therefore, 
all formulae describing the different types of the emission are 
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now completely different. The synchrotron and SSC emission 
scale in the same way 



^syn,2(v) = ^ syn ,l(v)| — 

r, 



Fssc,i(v) - Fssc,i(v) 



2.25 



2.25 



(32) 



this implies a decrease of the flux when F is decreasing. Also 
the EIC emission is decreasing with decreasing T 

F ex t,2(v) = F exttl (v)|^J . (33) 

Then, when the synchrotron cooling is dominant, all three ra- 
diative processes produce an observed flux which decreases 
with decreasing T. Note that, in this context, a decrease of T 
implies a comparison of two states of the jet, not a deceleration 
of a single source. 

In Fig. [2 we compare the emission of different sources, 
located at different distances from black hole that are travel- 
ing with different velocities. In this test we assumed the fol- 
lowing physical parameters for the fastest source: Ti = 20, 
R\ = 2.4 xlO 16 [cm],B\ = 7 [G],y u = 6xl0 2 , Q !( i = 8xl(T 5 
[cirT 3 s ]. We calculate the parameters for the other sources 
with T < Ti according to our prescriptions. Moreover, we 
assume that the sources are traveling within a constant radi- 
ation field, produced by a spherical shell (R &xt = 3 x 10 17 
[cm]). The emission of the shell is approximated as a black 
body with temperature r ext = 2 x 10 4 [K] and peak luminosity 
Vmax^ext^max) = [erg s -1 ]. We assumed a redshift z = 0.5, 
and the Hubble constant Hq = 72 [km/(s Mpc)]. Since we will 
later apply our model to low energy peaked blazars the cho- 
sen parameters are somehow representative of these class of 
objects. 

The reason why we parametrize the fastest source and not 
the slowest one is purely "technical". In order to calculate the 
particle distribution after the end of the injection we have to 
know the values of t/'„, U' and U'. The values of U'„ and 

ts syn exi ts 

U' ext can be calculated straightforwardly from our parametriza- 
tion, but the estimate of U' is more complex, since it depends 
on N(y), which in turn depends on U'. We then proceed by 
calculating the spectra for sources with different Lorentz fac- 
tors in a sequence, decreasing the value of T by a very small 
factor. This allows us to use the value of U' syD from the previ- 
ous step in the current computations if the change of U' syn from 
one step to another is small. This method requires that in the 
first step the synchrotron radiation field energy density is neg- 
ligible with respect to the other energy densities. Since (for our 
set of parameters) the fastest source cools mainly by the exter- 
nal radiation field, we always start the computations from the 
largest value of T. 

Our set of input parameters implies U'^ » U' B » U' 
for the fastest source (Fig.^d). Therefore, in agreement with 
our analytical predictions, initially Feic(v) remains constant 
whereas F syn (v) and Fssc(v) are increasing fast for decreas- 
ing F (Fig.^c)- However, the decrease of the Lorentz factor 
quickly reduces the importance of the external radiation field in 
the comoving frame. Thus, the synchrotron emission becomes 
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Fig. 1. In the upper panels we compare total emission of 
sources located at different distances from the center and trav- 
elling with different bulk Lorentz factors (F = 20 — » 11 panel 
a, and F = 1 1 — > 2 panel b). The total energy transferred to the 
electrons is exactly the same for each source. In this particular 
case the electrons are injected with a monoenergetic distribu- 
tion. Panel c) shows, separately, the fluxes of the different emis- 
sion mechanisms, integrated over narrow frequency ranges: 
synchrotron (10 n -10 12 Hz), SSC (10 14 -10 15 Hz) and the ex- 
ternal inverse-Compton emission (10 19 — 10 20 Hz). In the lower 
panel (d) we compare the different types of the energy densi- 
ties that control the particle cooling. This shows that, in this 
particular simulation, the external inverse-Compton scattering 
dominates the cooling process for large values of F whereas the 
synchrotron cooling is dominating for T < 10. 

the dominant cooling process (Fig.[2d). This leads to a fast de- 
crease of F eic(v) and also causes the decrease of F syn (v) and 
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Fig. 2. Top panels: comparison between spectra produced by 
different sources in the case of monoenergetic injection. In this 
particular test the SSC emission appears almost as important as 
the synchrotron radiation for sources with F ~ 10, for which the 
efficiency of the EIC process is significantly reduced. Note that 
also the second order SSC process appears significant in this 
particular case. The spectra of the different emission processes 
are presented in Fig.[3]for three representative sources. Bottom 
panels: as in Fig. 1, it is shown how the fluxes of the differ- 
ent emission processes, and the corresponding energy densities, 
scale as a function of T. 

^ssc(v), as predicted (Fig.^c)- Note that the parameters cho- 
sen for this particular simulation gives a very weak contribution 
of the SSC process to the total emission that appears visible 
only in the sources with very small T values (Fig^b). 

In the second test (Fig.|2} we modified some of the physical 
parameters (R[ = 2.4 x 10 16 -> 7.5 x 10 15 cm, B' = 7 -> 5 G, 
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Fig. 3. The components in the total spectrum produced by three 
different sources characterized by different values of T. The 
spectra are chosen from the sequence presented in Fig. |2 In 
the fastest source (F = 20) the external inverse-Compton emis- 
sion (EIC, long dashed line) appears dominant. In the moder- 
ately fast source (F = 11) the synchrotron (dotted line), SSC 
(short-dashed line) and the EIC processes are producing simi- 
lar amount of radiation. Moreover, also the second order SSC 
process (dash-dot line) is not negligible for this case. A rela- 
tively slow source (r = 2) is producing mostly synchrotron and 
SSC emission. 

Qi \ = 8 x 10~ 5 — > 1.5 x 10~ 2 crrT 3 s~') in order to amplify the 
SSC emission. However, even for this new set of the parameters 
the SSC and synchrotron cooling are negligible in the fastest 
source (Fi = 20) where the EIC emission dominates as in the 
previous test (Fig.|2]a). As we already mentioned the EIC cool- 
ing decreases quickly with the decrease of F at the advantage of 
other cooling processes. In this particular case, when F ~ 10, 
the synchrotron and SSC cooling processes become dominant 



and almost equally important ({/' 



^syn 



» c/; xt , Fig.EJd). 



Moreover, also the second order SSC emission appears signif- 
icant in this particular test. However, the energy density of the 
SSC emission (U' ssc ) that describes the particle cooling due to 
2 nd order SSC is for all the sources a few times smaller than U' B 
(Fig.|5]c). Note that if the SSC emission were stronger than the 
synchrotron radiation (U' > U' B ), the 2 nd order SSC scattering 
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would occur in the Klein-Nishina regime. This would require a 
special treatment of the cooling process and a numerical solu- 
tion of the kinetic equation. Our parameter choice is such that 
U' syn < U' B in all cases. As a consequence, the numerical calcu- 
lations are always in perfect agreement with the analytic pre- 
dictions. In Fig.|3]we show, separately, the spectra produced by 
the different emission mechanisms, and how their importance 
changes varying F. The fastest source (F = 20) is emitting the 
dominant part of the available energy through the EIC scatter- 
ing. The moderately fast source (F = 11) uses the synchrotron 
and the SSC processes to radiate the particle energy. The EIC 
in this source is significantly less important, producing a flux of 
the same order of the second order SSC emission. Finally, the 
slowest source (F = 2) is producing mostly synchrotron and 
SSC radiation. The parameters of the sources span the ranges: 
R' = 8 x 10 15 -48x 10 13 cm, B' = 5 — > 1.6 x 10 4 G, and 
# = 7.1xl0 3 ^2.3xl0 7 cm. 

In all cases we calculate the optical depth for the pair pro- 
duction process inside the source due to the interaction between 
y-rays and lower energy photons. This process appears negli- 
gible for all sources. 

We estimated also the energy density of the radiation pro- 
duced by the accretion disk {U' disk ) in the comoving frame ac- 
cording to the prescription of Ghisellini & Madau ( 1996 1. For 
this estimate, we assume that the total luminosity of the disk is 
ten times greater than the luminosity of the external radiation 
field, which is a good approximation if we assume that the ex- 
ternal radiation field is produced by the broad line region. The 
distance between the putative black hole and the source (D) 
is calculated assuming that the jet is conical, with an aperture 
angle of 0.1 rad (~5 degrees), giving D = R/Q.l. The disk en- 
ergy density starts to dominate over U' ext at relatively small dis- 
tances to the center, corresponding to small T-values (r < 5). 
Nevertheless, U' disk remains always at least two orders of mag- 
nitude weaker than the energy density of the magnetic field and 
the internal radiation fields (U' and U' ssc ). Therefore, the disk 
radiation field is negligible as a cooling agent of the particles 
inside the source. 

The above tests show that according to our prescription 
the fast moving sources radiate away the dissipated energy 
through the inverse-Compton scattering off the external pho- 
tons, whereas the slow moving sources may efficiently produce 
synchrotron and also SSC emission. The relative level between 
the SSC and the synchrotron emission depends on the phys- 
ical parameters of the source. What is important is that the 
same amount of the dissipated energy is radiated away, but 
with different emission processes. We can then observe com- 
pletely different spectra originating in the same jet, according 
to whether the dissipation takes place in its inner (small re- 
values) or outer (large T) part. For instance, the same cosmic 
source can display a state with a moderate optical synchrotron 
emission and a strong y-ray EIC radiation (at MeV-GeV en- 
ergies), or a strong optical emission and a strong SSC radi- 
ation (at keV energies) without a corresponding higher level 
of the MeV-GeV flux. Dramatic variability of a particular ob- 
ject would not correspond to dramatic variations on the jet effi- 
ciency (we have assumed that the blobs always carry the same 
energy), but would rather reflect the place, along the jet, where 



dissipation is taking place (if this location is in turn linked with 
the bulk Lorentz factor). 

Our model is rather "economic" for the jet, in the sense 
that the total energy radiated during different flares is approxi- 
mately the same, even if the variations in specific bands can be 
dramatic. The fact that a synchrotron (in the IR-optical band) 
and an X-ray flare may not be accompanied by a simultane- 
ous flare in the MeV-GeV band is the main prediction of our 
model, that can be easily verified by the coming new y-ray ob- 
servatories (GLAST, AGILE). These kind of observations may 
prove if the jet is varying its (likely kinetic, but possibly also 
magnetic) power or if it is instead working at a constant average 
efficiency. 

The monoenergetic injection used in our test provides a 
very simple description of the model and allows to describe 
the scaling of different radiative processes in an analytic way. 
However, such simple injection is corresponding to a particle 
distribution which is a power law, of slope n = 2, which can- 
not explain accurately most of the observations. Therefore, in 
the next step we discuss the injection of a power law particle 
distribution. 

2.4. Injection of a power law electron distribution 

In a homogeneous source, the spectrum of the observed emis- 
sion depends directly on the shape of the particle distribu- 
tion. Therefore, to precisely reproduce the observed spectra of 
blazars, we need a particle energy spectrum that is more com- 
plex than a single power law provided by the monoenergetic 
injection. Quite a large flexibility in spectral fitting is obtained 
if we assume a power law energy distribution of the injected 
particles 



fi(r) = Qi 1 ". for 7i < T < Tn 



(34) 



For j\ > 1 and n > 2 such injection with the simultaneous 
radiative cooling provides a particle distribution that can be ap- 
proximated by two power law functions 



My) 



Ky- 2 
K 2 y- (n+1) 



for y rain < y < n 
for ji < y < y max 



(35) 



where K 2 = Ky" 1 and the particle density (K) is directly re- 
lated to the injection rate (Qi) and the cooling (C) efficiency 



K = 



a 



C(n - l)y[ 



(36) 



The minimum energy y^ is calculated in the same way as 
for the monoenergetic injection (Eq.I15> whereas y max is a free 
parameter, whose value is not important if n > 2. The injection 
rate can be calculated from the simple relation 



Gi 



T'V'iy^t-y -")/(! -n) 



(37) 



where N e is the total number of injected particles. 

We will assume that sources with different T-values may 
have different slopes and different energy ranges of the injec- 
tion, but we still maintain our main assumption on the total 
energy of the injected electrons, which must be the same for all 
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sources (Eq.0. This gives the following relation between the 
total number of particles, average energy and F in two different 
sources 



N, 



e,2 



N, 



(38) 



Using the above description we can calculate the emission 
from different sources assuming different injection profiles but 
using the same amount of the injected energy. Even if the power 
law injection is more realistic than the monoenergetic injection, 
the results obtained are qualitatively the same for both cases. 

We will directly apply our simulations, done for a power 
law injection case, to the recent observations of the blazar 3C 
454.3, which we will use as an illustrative example. 



3. Application to 3C 454.3 

Recent observations of the blazar 3C 454.3 showed a huge flare 
in the optical and X-ray bands (Fuhrmann et al. |2006 Pian 
et al.HHH Villata et al. l2006l Giommi et al. 12006b . We will 
investigate if such activity is in agreement with the prediction 
of our model, in which sources with relatively small velocities 
(r < 10) can produce mostly synchrotron (IR — > UV) and SSC 
(UV — > hard X-ray) emission. Therefore in the next step we 
will try to reproduce this unusual activity. 

In Fig.|4]we collected the multifrequency observations of 
3C 454.3. Most of the data are historical observations made at 
different epochs, illustrating the low level of the optical and X- 
ray emission. The historical data from the radio frequencies up 
to the optical range are taken mostly from the NED database. 
However, in this range we are plotting also a few points from 
other sources (see Pian et al. 120061 and references therein). In 
the X-ray range we show the low/hard emission level observed 
by BeppoSAX in 2000 (Tavecchio et al.EOOS- Finally, in hard 
X-rays and y-rays we show the COMPTEL and EGRET mea- 
surements made in the years 1991-1994 (Zhang et al. 120051 
Hartman et al. I1999i . The high level of the optical and X-ray 
emission is shown by the data taken simultaneously by three 
experiments (UVOT, XRT, BAT) onboard the Swift satellite 
and by the quasi-simultaneous ground-based data taken by the 
REM telescope (Giommi et al. 2006). We will focus onto two 
states of the source: state 1 corresponds to the high level of 
y-ray emission, low/hard X-ray flux and low IR-optical flux, 
while state 2 corresponds to the recent optical and X-ray flare. 
Note that position of the synchrotron and the IC peak as well, 
is not well determined for the second state. 

We first try to reproduce state 1 . Assuming that the y-rays 
are produced by the EIC scattering inside a relatively fast mov- 
ing source, we obtained the best fit using the following param- 
eters: rj = 5i = 11, R\ = 1.1 x 10 16 cm, B\ = 4.8 G, y u = 90, 
y max = 5 x 10 3 , = 25.4 cirT 3 s~\ n x =2.1. For the ex- 
ternal radiation field we used Vmax£ext(vmax) = 10 45 erg s , 
r ext = 2 x 10 4 K, = 2.8 x 10 17 cm. The spectral index 
of the injected particles is relatively soft («i = 2.1), to explain 
the soft spectrum observed in the IR-U V range and also the 
relatively soft spectrum in the MeV-GeV range observed by 
EGRET. 



In order to reproduce state 2 we calculate a sequence of 
spectra decreasing the value ofF(r= 11 — > F = 6.25) by 
steps of Ar = 0.25. Moreover, in each step we increase the 
spectral index of the injected particle distribution by a factor 
An = 1/3. Also the minimum energy of the injected particles 
y; increases by a factor Ay; = 0.6yi,i whereas y max is assumed 
to be constant. According to our main assumption the changes 
of the injection spectrum are compensated by the change of 
the normalization (Eq. 1371 . to provide the same amount of the 
injected energy for all states. Note that in this case we main- 
tain the viewing angle fixed (6 = 5.2°) and therefore, in gen- 
eral, 5 + T. The best fit for state 2 is reached for T2 = 6.25, 
6 2 = 9.42, R' 2 = 3.7 x 10 15 cm, B' 2 = 35.7 G, y i>2 = 6 x 10 2 , 
Q i 2 = 3.4 x 10 13 cirT 3 s -1 , « 2 = 5.26. Note that while the 
set of parameters used for state 1 was chosen to best represent 
the data, the parameters for state 2 are given by our assumed 
parametrization, and are not free. Note also that we here appro- 
priately calculate the synchrotron self-absorbed spectra, that 
we neglected previously. 

Self-absorption truncates the synchrotron spectrum below 
the IR range. This effect appears frequently in the modelling of 
the high energy emission of blazars and is related to the size 
of the source that must be relatively small (10 15 ^ 16 [cm]) in or- 
der to explain the observed variability time scales. Such a small 
source must be also relatively dense (both in particles and mag- 
netic field) to explain the observed emission level and therefore 
it is optically thick at radio frequencies. The radio emission 
should be produced by more extended radio structures, well 
visible in the VLBI maps. Moreover, this particular blazar does 
not show variability of the radio flux clearly correlated with the 
optical, X-ray and y-ray activity (e.g. Villata et al. 2006 1. This 
suggests that the quiescent radio and the variable optical-to-y- 
ray emission are produced in different regions of the source. At 
any rate, the analysis of the radio emission produced probably 
by an inhomogeneous jet is out of the scope of this work. 

Note that our model reproducing state 2 is very similar to 
the spectrum resulting from the modelling (of the same state) 
made by Pian et al. (2006 1. In that paper, the emitting region 
was assumed to be outside the broad line region, in order to 
neglect the external radiation and the corresponding external 
Compton emission. In order to compensate the relatively low 
values of the magnetic field and the particle density, appropri- 
ate at such large distances from the center, it was necessary to 
adopt a large T, giving a Doppler factor 5=15. This difference 
in the T-factors gives a potential tool to discriminate between 
the model presented here (which has 6 = 9.4) and the model 
discussed in Pian et al. (2006 ). In fact, in our model, the appar- 
ent velocities of VLBI knots (after an optical and X-ray flare) 
should be smaller than in the Pian et al. model. Unfortunately, 
the predicted differences are not large, and probably not mea- 
surable with enough accuracy for this particular blazar. 

The presented simulation shows that our approach can 
well explain two different levels of the emission observed in 
3C 454.3. We stress that what we proposed is the most "eco- 
nomic" way to reproduce dramatically different levels of emis- 
sion, even if the jet is using always the same amount of energy. 

The strong prediction of this approach is that the level of the 
y-ray MeV-GeV emission should be lower (or not particularly 
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Fig. 4. The multifrequency observations of 3C 454.3 for two different emission states together with the results of our modelling. 
The low level of the optical and X-ray emission (state 1) is represented by data (filled circles) obtained at different epochs 
(references to the data points in the text) whereas the high synchrotron and X-ray state (state 2) is represented by the observations 
made by the Swift satellite and the REM telescope (Giommi et al. 2006 1. The left panel shows a sequence of spectra with different 
Lorentz factors: the model with the largest F reproduces state 1 and the model with the smallest F can explain state 2. In the right 
panels we show how the different emission processes contribute to the total spectrum. 



bright) when the level of the IR-to-X-ray emission is high. 
For the moment there is no evidence for such behaviour. We 
must wait for simultaneous optical, X-ray and especially y-ray 
observations of blazars to verify if our idea is correct. 

4. Summary and conclusions 

We have proposed a model to explain different levels of the 
emission observed from the IR range up to the y-ray band 
in a powerful blazars. The main assumptions of our model is 
that the emission zone of the jet is dissipating always the same 
amount of energy, but that the emission zone can be located 
in different parts of the jet, farther from the central engine if 
the bulk Lorentz factor is larger. We have discussed a partic- 
ular scenario where this naturally occurs, namely the internal 
shock scenario. While this scenario is not mandatory to de- 
velop our main idea, it is the only existing scenario where we 
can work out how all the relevant quantities (emission volume, 
magnetic field, particle density and so on) scale with the bulk 
Lorentz factor F. The obtained results are certainly promis- 
ing, since even with our simplifying assumption (the emission 
sites always dissipate the same amount of total energy) we 
were able to explain dramatically different states observed at 
different frequency bands. In particular, we were able to ex- 
plain two very different states of the blazar 3C 454.3. This is 
achieved by noting that when the bulk Lorentz factor F is small, 
the shell-shell collision in the internal shock scenario occur 
closer to the apex of the jet, therefore in a more compact re- 



gion, with a larger magnetic field. Since the external Compton 
radiation is dimmed in this case (being less boosted in the co- 
moving frame), the electron energy is radiated more through 
the synchrotron and the self-Compton processes, therefore not 
at MeV-GeV energies, but in the IR-optical (through syn- 
chrotron) and X-rays (through self-Compton). Lowering the 
F-factor there is therefore a "transfer of power" from higher to 
lower frequencies. 

If this idea is correct, there should exist sources flaring 
in the optical and X-ray bands, but not in y-rays, where the 
flux can even decrease. This prediction can be easily tested by 
the forthcoming y-ray satellites, AGILE and GLAST, together 
with simultaneous optical and X-ray observations: had MeV- 
GeV observations been available during the 2005 flare of 3C 
454.3, such a test could have already been done. 

The possibility that the jet can dissipate part of its bulk ki- 
netic energy at small distances from its apex, and at small dis- 
tances from the accretion disk, may at first sight contradict what 
claimed in Ghisellini & Madau (1996), where strong dissipa- 
tion close to disk was excluded. But in that case "dissipation" 
meant strong production of y-rays, dominating the radiative 
output. Here, instead, when the dissipation occurs "early" (i.e. 
close to the disk) the external Compton flux is unimportant, 
and the emission at high energies, made by the second order 
self-Compton emission, is never dominant. Therefore, even if 
a fraction of y-rays get absorbed in y-y collision producing 
pairs (likely with the X-ray emission from the disk corona), the 
spectrum is not strongly reprocessed, and the emission from the 
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newly created pairs is not important enough to "fill the valley" 
between the two broad peaks of the blazar SED. 

We do not know the relative fraction of time a jet spends 
producing blobs with large and small bulk Lorentz factors. 
Consider also that our current knowledge of blazars may be 
biased by the blazars observed (and detected) by EGRET. But, 
as mentioned in the introduction, EGRET detected only one 
fourth of the radio brightest blazars. This may suggest that a 
small value of the bulk Lorentz factor is the rule, not the excep- 
tion, since in this regime the synchrotron emission is favoured 
with respect to the external Compton emission. The factor ~20 
better sensitivity of GLAST with respect to EGRET will surely 
be crucial to solve this issue, and to shed light on these partic- 
ular aspects of powerful relativistic jet. 

Acknowledgements. This research has made use of the NASA/IPAC 
Extragalactic Database (NED) which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under contract with the 
National Aeronautics and Space Administration. We acknowledge the 
EC funding under contract HPRCN-CT-2002-00321 (ENIGMA net- 
work). 

References 

Begelman, M„ C, Rees, M, J., & Sikora, M. 1994, ApJ, 429, 57 
Fossati, G., Maraschi, L., Celotti, A., Comastri, A. & Ghisellini, G., 

1998, MNRAS, 299 
Fuhrmann, L., Cucchiara, A., Marchili, N., et al. 2006, A&A, 445, LI 
Giommi, P., Blustin, A. J., Capalbi, M., et al. 2006, astro-ph/0606319 
Ghisellini, G, & Madau, P., 1996, MNRAS, 280, 67 
Ghisellini, G., Celotti, A., Fossati, G, Maraschi, L. & Comastri, A., 

1998, MNRAS, 301,451 
Ghisellini, G, 1999, Astron. Nachrichten,320, 232 
Guetta, D., Ghisellini, G, Lazzati, D., & Celotti, A., 2004, A&A, 421, 

877 

Hartman, R.C., Bertsch, D.L., Bloom, S.D., et al. 1999, ApJS, 123, 79 
Krawczynski, H., Hughes, S., Horan, D., et al., 2004, ApJ, 601, 151 
Pian, E., Foschini, L., Beckmann, V., et al., 2006, A&A, 449, L21 
Kardashev, N.S., 1962, Soviet Astronomy-AJ, 6, 317 
Lazzati, D., Ghisellini, G„ & Celotti, A., 1999, MNRAS, 309, L13 
Meszaros, P., 2006, to appear in Reports on Progress in Physics 

lastro-ph/0605208! 
Rees, M.J., 1978, MNRAS, 184, 61 

Sikora, M., Begelman, M.C. & Rees, M.J., 1994, ApJ, 421, 153 
Spada, M., Ghisellini, G., Lazzati, D. & Celotti, A., 2001, MNRAS, 
325, 1559 

Tavecchio, F, Maraschi, L., Ghisellini, G., et al. 2002, ApJ, 575, 137 
Villata, M., Raiteri, C. M., Balonek, T, J., et al., 2006, A&A, 453, 817 
Zhang, S., Collmar, W., & Schnf elder, V., 2005, A&A, 444, 767 



